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ABSTRACT 

A model to calculate fatigue life is developed based on the assumption that fatigue life is 
entirely composed of crack growth from an initial microstructural inhomogeneity. Specifically, 
growth is considered to start from either an ellipsoidal void, a cracked particle, or a debonded 
particle. The capability for predicting fatigue life from material microstructure is based on linear 
elastic fracture mechanics principles, the sizes of the crack-initiating microstructural 
inhomogeneities, and an initiation parameter that is proportional to the cyclic plastic zone size. 

A key aspect of this modeling approach is that it is linked with a general purpose probability 
program to analyze the effect of the distribution of controlling microstructural features within the 
material. This enables prediction of fatigue stress versus life curves for various specimen 
geometries using distributional statistics obtained from characterizations of the microstructure. 
Results are compared to experimental fatigue data from an aluminum alloy. 


INTRODUCTION 


Classically, the fatigue life of a metallic specimen or structure has been thought of in three 
stages: crack initiation, crack growth and final fracture. As crack measurement techniques have 
improved and researchers have examined small crack growth, microcracks have been 
discovered which grow from very early in a structure's life. This period of small crack growth 
can be 50 to 90% of the fatigue life. This suggests that fatigue life can be considered as a two 
stage process of growth and fracture. Due to the length of time that a structure spends in this 
early phase of fatigue crack growth it is important to accurately model this behavior. 

In metallic structures there is generally a hierarchy of features which can control fatigue. 
Various researchers have identified some of these features. In 2024-T4 and 7075-T6, Erhardt 
et. al. [1], indicate cracks started at surface inclusions; Wang [2] in a study of 2024-T3 
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identified inclusions and tool marks as key initiators; Lukasak and Koss [3] in a study of a 
particle reinforced MMC identified two failure modes associated with the reinforcing particle; 
Lankford and Kusenburger [4] linked high strength steel crack initiation to cavitation caused by 
particle debonding, ; Goto, et al. [5] identified particles and surface scratches as fatigue 
initiators in heat-treated 0.45% C steel; Fischmeister, et al. [6] determined pores and particles to 
be the initiation sites of a Ni superalloy; Magnusen, et al. [7] indicated that in 7050-T7451 pores 
and particles act as initiators. Fracture mechanics and damage tolerant methodologies have been 
shown to work well for rogue inspectable flaws on the order of 0.01-0.05 inches. The 
motivating force behind this work is to model the role of alloy microstructure on fatigue 
durability as identified by cracking that originates from a feature of the microstructure 
approximately 0.0005-0.01 inches in diameter. The crack is then grown to an inspectable 
dimension of economic consequence (i.e., a crack requiring diagnostic or corrective action 
during the operational life of a part). 

The principal types of microstructural features of concern are micropores and constituent 
particles. These microstructural inhomogeneities act as sites for the initiation of fatigue cracks. 
The next section discusses the basic fracture mechanics approach used to model early fatigue 
growth. Following sections then discuss the infinite life stress sensitivity and sensitivities of 
the model to material parameters. Finally, life predictions using the model are compared to 
experimental lives on an aluminum alloy. 


MICROSTRUCTURAL BASED FATIGUE MODELING 


The method used for the fatigue life analysis is based on the assumption that fatigue life is 
entirely composed of crack growth from an initial inhomogeneity [8]. Specifically, growth is 
considered to start from either an ellipsoidal void, a cracked particle, or a debonded particle. It 
is assumed that each type of inhomogeneity contains a pseudo crack of equatorial length b. As 
shown in Figure 1, the particle or void height is h and its width is 2R. The remote stress, a, is 
applied in the z direction, perpendicular to the plane of crack growth. It is assumed that for a 
given far-field cyclic stress, Aa, the pseudo crack forms in the first few cycles. The pseudo 
crack can be considered a perturbation of the surrounding stress field for the purpose of 
calculating the available energy for crack extension in the vicinity of a local pore/particle stress 
concentration. The length of this pseudo crack is assumed to be proportional to the cyclic plastic 
zone size ahead of the crack as illustrated in Figure 1. To make these ideas more precise, 
consider the Trantina-Barishpolsky effective stress intensity factor [9], for the flaw geometry 



shown. The stress intensity factor associated with the crack surrounding the inhomogeneity can 
be expressed as. 


Kt— B — PoV^tb 

where the dimensionless term [3 addresses the local microfeature and is given by. 
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Here B is a constant that has the value 1 for a void, 2 for a bonded cracked particle, and 0.3 
for an debonded particle; k t is the local elastic stress concentration factor for the ellipsoidal void 
or particle without the crack. Note that p depends on the geometry through the ratio R/(b+R), 
so that for b » R, Equations 1 and 2 reduce to the stress intensity factor solution for an 
embedded elliptical crack, where the void or particle lacks influence on the stress intensity 
factor. At the other extreme, for b = 0, from Equation 1, the stress intensity factor is zero. 
Under the imposition of a far-field cyclic stress range Act, the cyclic stress intensity factor is, 

AK T — b = PAaVib (3) 

An estimate of the cyclic plastic zone size associated with AKt-b is, [10] 
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where a y is the yield strength of the material. The procedure used to calculate the parameter, b, 
is to equate it to this measure of plastic zone size. Thus, solve Equation 3 for b and equate it to 
the plastic zone size of Equation 4. This amounts to finding the ratio R/(b+R) that satisfies 
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THRESHOLD STRESS FOR INHOMOGENEITY INITIATED FATIGUE 


Equation 5 can be used directly to estimate the fatigue threshold stress, or the stress which 
yields infinite life, AOoo- For infinite life, b = 0, and 
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Although the number of cycles to failure is sensitive to the initial void or particle size, 
Equation 6 suggests that if fatigue life is controlled by the growth of a crack from an 
inhomogeneity, then the infinite life stress is independent of void or particle size. It is 
proportional to yield strength and shows a strong dependence on local stress concentration 
factor, kt 5 and the type of inhomogeneity through B. For voids, B — 1 and Equation 6 
becomes. 
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For cracked particles, B = 2 and Equation 6 becomes, 
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A plot of the sensitivity of the infinite life stress to the local stress concentration factor is 
shown in Figure 2 for both voids and cracked particles. Figure 2 shows that voids affect the 
infinite life stress more severely than uncracked particles. For voids, a lower limit on the stress 
concentration factor is a smooth sphere for which k t is about 2. At this point the ratio of infinite 
life stress to yield stress is about 0.9. For Equations 7 and 8 it is assumed that the yield 
strength is an upper limit on the infinite life stress for the mechanism of fatigue crack growth 
from a void or particle. This assumption sets limits on the stress concentration factors below 
which little improvement in infinite life stress is expected. For voids, improvements are 
possible up to the point at which the void is the shape of a smooth sphere which is the minimum 
stress concentration possible. For non-spheres there is a dramatic drop in the infinite life stress 
for small increases in the local stress concentration. For a stress concentration of 3, the ratio of 
infinite life stress to yield strength has dropped to about 0.6. This corresponds to a circular 
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cavity of elliptical cross-section with aspect ratio 2.5, as depicted in Figure 3. For particles, the 
value of stress concentration factor at which the infinite life stress begins to decrease is 1.62. 
The stress concentration of the particle depends both on particle shape and particle elastic 
modulus. For a bonded uncracked particle, stress concentrations in the matrix are typically less 
than 2 [1 1]. A plot of the stress concentration factor for an isolated rigid elliptical cylindrical 
particle is shown in Figure 4. For far-field uniform tension, the stress concentration factor for 
tangential tension is less than 2 for aspect ratios ranging up to 30. This suggests that particle 
cracking will typically be required before particles can control the infinite life stress. In 
previous studies on the effects of cracked particles [12,13], it has been shown that cracked 
particles behave like high aspect ratio voids. This is expected to severely limit the infinite life 
stress, as shown in Figure 2. 

Since the stress concentration factor, k^ is taken to be a local value at the point where a 
crack starts to grow. Equation 6 shows that the infinite life stress is sensitive to the smoothness 
of the void or particle. 


SENSITIVITY OF FATIGUE LIFE CURVE TO MODEL PARAMETERS 


Crack growth emanating from a metallurgical inhomogeneity is assumed to follow a power 
law, 


da 
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(9) 


where ^ is the crack growth rate, AK is the cyclic stress intensity, and A and m are 

material constants. For the aluminum plate used in this study, the coefficients obtained by 
fitting the equation to long crack growth data were, A = 3.9X10" 9 , m = 4.18 for an R-ratio of 
0.1. AK is in units of ksi ^[m. 

The sensitivity of the fatigue model to changing material and microstructural parameters can 
be assessed by calculating the cyclic stress versus life (S/N) curves for the varying material 
conditions. Specimen failure occurs when K max = Ki c ; a value of Kj c = 25 ksi Vin was used. 
The influence of particle size on fatigue life is shown in Figure 5 for particle sizes ranging from 
0.001 in. to 0.01 in. at a stress concentration value of 2. The yield strength was fixed at 70 ksi. 
The effect of increasing particle size is to decrease life throughout the range of stress levels. 
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The infinite life stress, however, is insensitive to particle size. Similar trends hold for 
increasing void size at constant aspect ratio. The sensitivity of fatigue life to yield strength is 
shown in Figure 6. Here the particle size is fixed at 0.005 in. and the yield strength is varied 
from 60 ksi to 80 ksi. As seen, there is almost no change in fatigue life with stress level above 
the infinite life stress, AOoo. The effect of yield strength on AOoo is strong; for this mechanism, 
increasing yield strength increases the threshold stress for fatigue. The results shown in 
Figures 5 and 6 show that if simultaneously particle (void) size is decreased while increasing 
yield strength, fatigue life is expected to increase at all stress levels. The inhomogeneity aspect 
ratio or stress concentration has a similar effect on the fatigue life curve as yield stress does. 
Figure 7 shows a comparison of fatigue S/N curves for particles and voids at different shapes 
for a yield strength of 70 ksi. The particle or void size is fixed at 0.005 in. It is seen that 
increasing the stress concentration has a large effect on decreasing the infinite life stress. 


COMPARISON OF MODEL CALCULATIONS TO TEST DATA 


Open hole fatigue tests were performed to quantify the influence of the various 
microstructural features on fatigue performance of a high strength aluminum alloy. The fracture 
surfaces from the failed specimens were then examined in a scanning electron microscope and 
the size, location and type of crack initiating inhomogeneity were recorded. Failures in material 
A were controlled by microporosity, while the failures for material B were generally controlled 
by cracked particles located at or near the comers of the bore hole. Despite the change in 
initiation mechanism from material A, where micropores dominated, open hole lifetimes of the 
two materials are similar. 

The model described in the previous sections has been used to calculate open hole fatigue 
performance of these materials. These calculations have been made using data obtained from 
the failed fatigue specimen fracture surfaces. Additional input for the calculation included the 
yield strength, crack growth properties of the material and the local stress concentration factor at 
the microstructural feature. The local stress concentration factor was determined empirically 
from fitting the model with the infinite fatigue life data for each material variant. The values 
obtained for the local stress concentration factor for material A and B are 7.5 and 4.25, 
respectively. These values can be interpreted to say that the local stress concentration factor 
became less severe as the microstructural feature changed. The growth model was linked to the 
program PROBAN [14] to permit use of statistical information on particle and pore 



distributions. PROBAN is a general purpose probability analysis program which can be linked 
to user-defined subroutines to predict crack growth. 

The calculated stress versus life curves for each material variant, along with the open hole 
fatigue test data are shown in Figures 8 and 9. The model was used to calculate the 5%, 50%, 
and 95% fractiles for each material. (The curves shown are not smooth because only a limited 
number of stress levels have been analyzed with the model.) In general, the results of the model 
are in good agreement with the test data. Additional refinement of the model is expected to yield 
improvement in the ability to predict fatigue behavior of the material variants. 


CONCLUSIONS 


• Fracture mechanics principles can be applied to microstructural features (0.0005-0.01 inch) 
that govern fatigue life for both smooth axial and open hole specimen geometries. 

• The effects of the initial inhomogeneity size, shape, and location distributions on specimen 
fatigue life distributions can be analytically predicted. 

• Calculations of specimen fatigue S/N curves agree well with actual test data. 

• Typically cracked particles are required before particles will control the infinite life stress. 

• A particle’s or pore's local stress concentration factor is important. The infinite life stress is 
a maximum when the local stress concentration factor is small. 


FUTURE WORK 


A goal of modeling is to predict material performance based solely on characterization of the 
micro structure rather than a postmortem analysis of fatigue failures. Therefore, it is necessary 
to take metallographic characterizations of the microstructure obtained on random planes and 
scale the distributions of controlling features to obtain the extreme distributions that are 
observed through fractography. Extreme value statistics will be used to establish the 
relationship between the random plane characterizations and the extreme value distributions of 
constituent particles or voids. Future work will develop capability to predict material fatigue 
performance based on material microstructure. 
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Figure 2. Sensitivity of threshold stress to stress concentration. 
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Figure 3. Stress concentration for a circular cavity of elliptical cross-section. 
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Figure 4. Stress concentration factors for isolated rigid, elliptical, cylindrical inclusion 
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Figure 5. Sensitivity of fatigue lifetime curve to particle size for yield strength of 70 ksi. 
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Figure 6. Sensitivity of fatigue lifetime curve to yield stress for particle size of 0.005 in., k = 2. 
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Figure 7. Comparison of fatigue lifetime curves for particles and voids at different particle shapes 
for yield strength 70 ksi., particle size = 0.005 in. 
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8. Calculated stress versus life curves for Material A along with the open hole fatigue 
test data. Tests conducted at R=0.1, LT orientation, freq =30 Hz, lab air. 
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Figure 9. Calculated stress versus life curves for Material B along with the open hole fatigue test 
data. Tests conducted at R=0.1, LT orientation, freq.=30 Hz, lab air. 
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